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CONSTRAINTS ON THE EMISSION AND VIEWING GEOMETRY OF 
THE TRANSIENT ANOMALOUS X-RAY PULSAR XTE J1810-197 



OO 

O 

o 

(N 



43 

6 



> 

O 

(N 

cn 
o 

oo 
o 



- 1—1 

x 



Draft version March 13, 2008 

ABSTRACT 

The temporal decay of the flux components of Transient Anomalous X-ray Pulsar XTE J18 10-197 following 
its 2002 outburst presents a unique opportunity to probe the emission geometry of a magnetar. Toward this goal, 
we model the magnitude of the pulsar's modulation in narrow spectral bands over time. Following previous 
work, we assume that the post-outburst flux is produced in two distinct thermal components arising from a 
hot spot and a warm concentric ring. We include general relativistic effects on the blackbody spectra due to 
gravitational redshift and light bending near the stellar surface, which strongly depend on radius. This affects 
the model fits for the temperature and size of the emission regions. For the hot spot, the observed temporal and 
energy-dependent pulse modulation is found to require an anisotropic, pencil -beamed radiation pattern. We are 
able to constrain an allowed range for the angles that the line-of-sight {ip) and the hot spot pole (£) make with 
respect to the spin-axis. Within errors, this is defined by the locus of points in the ^-^-plane that lie along the 
line (£ + (3(R))(ip + /3(R)) »s constant, where (3(R) is a function of the radius R of the star. For a canonical value 
of R = 12 km, the viewing parameters range from ip = £ = 37° to (ip,Q = (85°, 15°). We discuss our results in 
the context of magnetar emission models. 

Subject headings: pulsars: individual (XTE J18 10-197) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) are peculiar high-energy 
pulsars whose observed luminosity greatly exceeds that which 
can be supplied by their rotational energy losses. These 
pulsars occupy a narrow range of spin periods (P ~ 2— 12 
sec) and are spinning down rapidly compared to the rotation- 
powered pulsars. For the vacuum dipole model, the timing 
properties of these pulsars imply an enormous magnetic field, 
B ~ 10 14 - 10 15 G. These relatively rare objects (~ 10 com- 
pared to ~ 1700 catalogued radio pulsars), generally display 
sinusoidal modulation in their pulsed flux, with a wide range 
of amplitudes (~ 10-80%) and are likely young (< 10 4 yrs), 
as more than half are associated with supernova remnants (see 
Kaspi 2007 for a recent review). AXPs can be understood 
within the context of the magnetar model developed by Dun- 
can & Thompson (1992) to explain the burst phenomenology 
of Soft 7-ray Repeaters (SGRs). The excess emission from 
both AXPs and SGRs, collectively referred to as magnetars, is 
powered by the decay of their extreme magnetic fields. This is 
suggested by the relatively high temperatures of their thermal 
emission (kT w 0.4-0.7 keV for blackbody fits), and frequent 
rapid (< 0.1 s) burst activity. The geometry and the proper- 
ties of the observed emission from the magnetars is of great 
interest for understanding how this activity arises. 

The recent discovery of an AXP fading from a long duration 
outburst offers the unique opportunity to probe the magnetar 
emission geometry evolution during this event. The P = 5.54 s 
Transient AXP (TAXP) XTEJ1810-197 was discovered in 
January 2003 by Ibrahim et al. (2004) using the Rossi X- 
ray timing explorer (RXTE) following a large eruption. Sub- 
sequently its flux decayed exponentially (t w 900 d) nearly 
back to a quiescent flux level as determined from serendip- 
itously archival X-ray observations (Gotthelf et al. 2004). 
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The earlier measurements indicate that the TAXPs 2-10 keV 
flux had increased by two orders of magnitude. However, of 
great interest, the quiescent luminosity is 100 times lower than 
for the well-established AXPs and SGRs, suggesting a large 
unidentified population of neutron stars (Gotthelf et al. 2004). 
In contrast, the magnetic field strength of XTE J18 10-197, 
B = 3 x 10 14 G as inferred from its spin-down properties, is 
typical of the magnetars. 

The flux and pulse evolution of XTE J 1 8 1 0-197 were mon- 
itored with the XMM-Newton X-ray observatory at roughly bi- 
yearly intervals starting Sept. 2003, yielding a total of seven 
epochs through Mar. 2006. The complete set of observations, 
together with their spectral modeling and interpretation, is 
described in detail by Gotthelf & Halpern (2007), with the 
earlier observations reported in Gotthelf & Halpern (2005) 
and Halpern & Gotthelf (2005). While analysis of phase- 
averaged spectra alone cannot distinguish among competing 
models for the AXP emission type and geometry, the addition 
of the steady change of the spectrum and pulse modulation 
over time greatly increases the diagnostic power. 

In this paper we present a detailed model for the en- 
ergy dependent pulse phase from XTE J18 10-197. This 
model accounts for the viewing geometry and surface emis- 
sion distribution. We include the general relativistic effects 
of light deflection and gravitational redshift and allow for 
anisotropic emission. We apply this model to a set of X-ray 
data acquired during the temporal evolution of the flux from 
XTE J18 10-197. This allows us to constrain the underlying 
emission geometry and radiation properties of this transient 
magnetar. 

2. TIME-DEPENDENT FLUX MODELING OF XTE J1810-197 

2.1. Model Motivations 

Since their discovery, spectra of magnetars have been fit- 
ted with a variety of models, generally including two com- 
ponents, such as blackbody plus power-law, atmosphere plus 
power-law (Perna et al. 2001a; Skinner et al. 2006), or ther- 
mal plus resonant cyclotron scattering (Rea et al. 2007), or 
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more recently, a magnetized atmosphere model with the in- 
clusion of scattering (Guver et al. 2007). From the point of 
view of phase-averaged spectral analysis alone, these models 
are generally statistically acceptable, and therefore none can 
be ruled out a priori. 

Phase-resolved modeling of the observed modulation, how- 
ever, can provide a much stronger constraint. This is partic- 
ularly true for the case of XTE J18 10-197 , given the wealth 
of data available at different epochs, while the object is cool- 
ing. From the point of view of being able to reproduce the 
observed time-dependent energy behavior of the pulsed frac- 
tions, the double blackbody model (made up of a hot and a 
warm component) put forward by Gotthelf & Halpern (2005) 
is quite promising, at least from a qualitative point of view. 
In fact, if, in the same energy band, the warm component is 
less pulsed than the hot one, then, as the hot component drops 
in flux faster than the warm one (as found in their fitting), 
the pulse fraction will tend to decline, as observed. Also, in 
their modeling, the area of the warm component, which be- 
comes more luminous than the hot component after the fourth 
epoch, increases at later times. 3 This, again, tends to produce 
a decline in the pulse modulation and counteract the increase 
that would otherwise have, due to the decreasing temperature. 

While it is tempting to model the thermal components with 
detailed magnetized atmospheres (e.g. van Adelsberg & Lai 
2006), these models might be problematic for the case of 
XTE J18 10-197 following its outburst. Invoking dissipation 
of a twist in the magnetic field lines (Beloborodov & Thomp- 
son 2007), the field in the emitting region is likely to have 
significant non-normal components. While the "twist model" 
nicely predicts the timescale of the outburst decay, the non- 
normal surface magnetic fields ( i.e. magnetic fields that 
emerge from the star surface at oblique angles) have yet to be 
fully realized in the magnetized atmosphere models. Since the 
predicted amplitude of the flux modulation strongly depends 
on the local magnetic field direction, by assuming a magne- 
tized atmosphere model for the thermal components we would 
introduce an a priori bias in our results. 

In this work, we prefer to take a more empirical approach by 
starting with the distribution of the emitted radiation over the 
stellar surface and allowing a degree of anisotropy (beaming 
factor) in the thermal (blackbody) components (following the 
methods of Pechenick et al. 1993; DeDeo et al. 2000; Perna 
et al. 2001b), and leaving the beaming factor as one of the 
model parameters. This approach turns out to be very valu- 
able in that we can use the energy-dependency of the pulsed 
fractions, together with their variation with time as the emit- 
ting region cools, to constrain both the viewing geometry and 
the beaming properties of the radiation simultaneously (§3). 
Our results can therefore be used as a guide for further theo- 
retical modeling aimed at understanding the detailed mecha- 
nisms that produce the observed thermal radiation. 

2.2. The Emission Model 

Our method for parameterizing the surface emission from 
XTE J 1 8 1 0-197 follows the example developed by Pechenick 
et al. (1993) with some generalizations. For our numerical 

3 This area increase might either be real, or a result of the fact that the 
warm component cannot be straightforwardly separated in the fits from the 
underlying surface emission, to which it approaches at later times. Either 
way, the pulse modulation in this model is expected to decrease with time, 
since the surface emission of the neutron star, if its temperature distribution 
traces that of a dipolar magnetic field, is not expected to be highly modulated 
(DeDeo et al. 2001). 




FIG. 1 . — Emission geometry on the surface of the neutron star (NS) for 
the model presented herein. A hot spot of temperature 7), and angular size 
Pi, is surrounded by a warm ring of temperature T w with outer angular size 
(3 W . As the neutron star (NS) rotates with angular velocity n(f), the angle 
a(t) is a function of the phase angle 7(?) = Q(t)t, and the angles i/> and £ 
between spin axis and viewing vector, and between spin axis and emission 
pole, respectively. Due to general relativistic effects, a photon emitted at a 
colatitude 9 on the star's surface which reaches the observer must be emitted 
at an angle S with respect to the star's normal at that point. 

modeling we consider emission from a hot spot of tempera- 
ture Th and angular radius surrounded by a warm ring of 
temperature T w and outer radius /3 W 

We indicate with a the angle that the axis of the hot spot 
makes with respect to the line of sight. This depends on the 
phase angle 7(f) = ft(t)t, as the star rotates with angular veloc- 
ity ft(t). If £ is the angle between the spot axis and the rotation 
axis, and ip the angle between the observer's direction and the 
rotation axis, then the angle a is given by 

a(t ) = arccos(cos ip cos £ + sin tp sin £ cos 7(f))- ( 1 ) 

The geometry is illustrated in Figure 1 . The surface of the 
star is described by the angular spherical coordinates (6,4>), 
and the coordinate system is chosen so that the z axis is in the 
direction of the line of sight to the observer. This is a nat- 
ural choice of coordinate system for our problem, since the 
observed flux is produced by all the photons that reach the 
observer at 00 along that axis (e.g. Pechenick et al. 1993; 
Page 1995). Note that, if there were no general relativistic ef- 
fects, a photon emitted at colatitude 9 on the star would only 
reach the observer if it were emitted at an angle 5 = 9 with the 
normal to the surface of the star. Because of general relativis- 
tic effects, however, a photon emitted at a colatitude 9 will get 
to the observer only if emitted at an angle 8 with respect to 
the surface normal (see Fig.l), where the relation between the 
two angles is given by the ray-tracing function 4 (Pechenick et 
al. 1993; Page 1995) 

9(S) = l RJ2R xdu I (|) 2 -d-2«)«V , 

(2) 

4 To improve the computational efficiency of the above equation we use 
the approximation presented in Beloborodov (2002). 



3 



having defined x = sin 5. Here, R/R s is the ratio of the 
stellar to Schwarzschild radius, R s = 2GM/c 2 (we assume 
M= 1.4M ). 

The hot spot is described by the conditions: 



<Ph, 



if a = 



and 



where 



a~Ph < 9 < a + (3 h 

2TT-(j) h p if a ¥0 and (3 h < a 



■■ arccos 



cos (3h - cos a cos 9 



sin a sin f 



(3) 



(4) 



(5) 



On the other hand, it is identified through the condition 

9 <9^(a,f3 h ,4>), if a=/0 and f3 h > a , (6) 

where the outer boundary 6^(a, /?/,, (f>) of the spot is computed 
by numerical solution of the equation 



cos /?/, = sin 61 sin a cos <\> + cos 6* cos a 



(7) 



Similarly, the warm ring is described on the star surface 
through the conditions 



Ph<0< /?„,, if a = 



(8) 



and 



a-0 w < 9 < a+0 w and a + /3 h <9 <a-(3 h 
* < d)<2n-^' 



2tt- 



< , 



and 



if a and f3 w < a. In the above equation, (f> p has the same 
functional form as </>£ in Eq (0, except for the substitution 
Ph — * fiw Finally, if a ^0 and f3 w > a, the ring is identified 
by the condition 



(10) 



where, again, the outer boundary of the ring 9™{a,f3 w ,<fr) is 
found by numerical solution of Eq. (|7), but with the replace- 
ment j3 h -> P w . 

In the following we assume a blackbody emission model 
for both the hot spot and warm ring, characterized by a uni- 
form temperature (7), or T H ,, respectively) over their stellar sur- 
face. As discussed above, we allow for the radiation from the 
two regions to be anisotropic, and parameterize the beaming 
of their local emission through the functions fh(5) oc cos"''(<5) 
(hot spot) and f w (S) oc cos"' 1 (5) (warm ring). This choice was 
initially motivated by the consideration that the hot spots are 
likely associated with regions of larger conductivity, where 
the magnetic field lines would be close to perpendicular to the 
surface of the star. This would produce an enhanced emissiv- 
ity at small 5. Our analysis (§3.2) then confirmed the validity 
of this choice by demonstrating that the modulation of the 
hottest region does indeed require a pencil-type anisotropic 
beaming pattern. 

The observed spectrum as a function of phase angle 7 is 
then obtained by integrating the local emission over the ob- 
servable surface of the star, accounting for the gravitational 
redshift of the radiation (Page 1995) 



2tt fl 2 



2xdx 



2tt 



d([> 
2^ 



k(6,&n[E 00 e- A °;T(6,4,)] 



(11) 



in units of photons cm" 2 s" 1 keV" 1 . In the above equation, the 
radius and energy as observed at infinity are given by R^ = 
Re~ As , and = Ee A ', where R is the star radius, E is the 
energy emitted at the star surface, and we have defined A s as, 



R 



(12) 



For the spectral function, given by n(E,T) = 
l/[exp(E/kT) - 1], the temperature T(6,(f)) is equal to 
Ti, if { 9, </>} satisfy any of the conditions ([3]) through (0, 
and it is given by T(9,<p) = T w if any of the conditions (HJ 
through ( [Tol l holds true. Correspondingly, the weighted 
intensity Iq{9, </>) is given by the functions //, or f K depending 
on whether the variables {9,4>} are inside the hot or warm 
region, respectively. 

The phase-averaged flux is readily computed as -F a ve(£oo) = 
1/27T / 2lr djFiEocj). The phase dependence 7 in Eq.dTTl 
comes from the viewing angles implicit in a(t ) and from the 
series of conditions d/3j through ( TTOb . Note that, as the star 
rotates, the only angle on which the flux depends is a(t), the 
angle that the line of sight makes with the axis of the spots. 
Since, in the magnetar model, the spots are likely to be corre- 
lated to regions with an enhanced magnetic activity, the angle 
a(t) can also be considered as the (phase-dependent) angle 
between the line of sight and a magnetically active region on 
the star during and following the outburst. When the star re- 
turns to quiescence, the temperature distribution on the star 
will reflect the overall magnetic field configuration. For most 
AXPs, the quiescent emission cannot be produced by a tem- 
perature distribution following a dipolar magnetic field (De 
Deo, Psaltis & Narayan 2000); in the case of XTE J1810- 
197, a detailed study of the quiescent emission, once the con- 
tamination from the heated regions has completely subsided, 
will be able to determine the detailed structure of the surface 
temperature distribution, and hence reconstruct the magnetic 
field structure in quiescence. This study will be performed 
in a forthcoming paper. However, a preliminary investiga- 
tion of the softest energy band in the latest data set (where 
the surface emission from the rest of the star is likely to be 
dominant) shows that the maximum of the pulsed emission 
remains in phase with the maximum in the hardest energy 
band (still dominated by the heated region). This result shows 
that the maximum of the quiescent emission comes from the 
region where the outburst occurred. Therefore, the active re- 
gion is likely to be associated with an enhanced magnetic field 
strength in quiescence. If the underlying B field is dipolar (or 
close to such), then the spot axis in our paper also represents 
the dipole magnetic field axis, and therefore the angle a(t) 
with respect to the z-axis (observer's viewing direction) would 
naturally be associated with the angle that the magnetic axis 
forms with respect to our line of sight as the star rotates. 

3. MODELING THE MULTI-EPOCH MODULATION 

The above model was coded in FORTRAN and fully im- 
plemented as an additive model with 10 parameters in the 
XSPEC spectral fitting software (Arnaud 1996). Specifically, 
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TABLE 1 

Spectral Results as a Function of the NS Radius - Model Fits for 
N H = 6.8 x 10 21 cm 2 , D = 3.3 kpc, i/> = £ = 7 = 0° (see text for 
definitions) 



Parameter 


2003 


2004 




2005 


2006 




Sep 8 


Oct 12 


Mar 11 Sep 18 


Mar 18 


Sep 20 


Mar 12 








# = 9 km; x 2 (dof) 


= 1.09(1914) 






l-T, (hf*\/\ 

Kl h (Kevj 


n qi 


0.95 


0.93 


0.91 


0.84 


n 71 


U . OU 


kT w (keV) 


0.32 


0.38 


0.35 


0.35 


0.31 


0.26 


0.23 


A (deg) 


10.06 


8.47 


6.69 


5.33 


3.81 


3.52 


4.08 


A< (deg) 


56.3 


38.3 


39.3 


42.3 


42.3 


55.1 


82.1 






R= 10km; x 2 (dof) 


= 1.11(1914) 






/-X. ftpAA 

kj;, (Kevj 


n qi 


0.95 


0.890 


0.91 


0.76 


U.D / 




fcr„, (keV) 


0.36 


0.40 


0.33 


0.35 


0.27 


0.24 


0.22 


A (deg) 


8.21 


7.21 


5.83 


4.46 


4.03 


3.12 


3.11 


Pw (deg) 


OC Q 
JJ.O 


29.4 


34.6 


30.2 


43.8 


j i .y 


DO. J 






R= 11km; xJ(dof) 


= 1.15(1914) 






Klh (KeV) 


n so 

U.07 


0.90 


0.87 


0.83 


0.75 


U.D / 


U.jO 


kT w (keV) 


0.35 


0.37 


0.32 


0.29 


0.27 


0.24 


0.21 


A (deg) 


7.05 


6.98 


5.25 


4.37 


3.44 


2.85 


3.27 


/Q (J^„\ 

Av (deg) 


31.3 


27.8 


31.5 


36.3 


37.1 


43.8 


61.3 






R = 12 km; x 2 (dof) 


= 1.09(1914) 






/t7], (keV) 


0.85 


0.86 


0.86 


0.82 


0.71 


0.63 


0.59 


W w (keV) 


0.32 


0.34 


0.32 


0.29 


0.25 


0.24 


0.22 


A (deg) 


7.06 


6.39 


5.01 


4.22 


3.45 


2.53 


2.27 


Av (deg) 


33.8 


28.1 


28.6 


32.1 


37.8 


41.2 


45.3 






#=13 km; xJ(dof) 


= 1.10(1914) 






kT h (keV) 


0.83 


0.85 


0.82 


0.76 


0.69 


0.61 


0.53 


W„, (keV) 


0.31 


0.33 


0.29 


0.26 


0.24 


0.23 


0.20 


A (deg) 


6.35 


5.84 


4.81 


4.14 


3.29 


2.44 


2.83 


Av (deg) 


31.4 


26.5 


29.6 


36.7 


37.9 


39.2 


49.2 






#=14 km; x£(dof) 


= 1.14(1914) 






feT* (keV) 


0.80 


0.81 


0.79 


0.74 


0.67 


0.59 


0.52 


kT w (keV) 


0.28 


0.30 


0.27 


0.23 


0.23 


0.22 


0.19 


A (deg) 


6.16 


5.79 


4.49 


4.07 


3.23 


2.40 


2.45 


Av (deg) 


33.26 


28.9 


30.7 


43.5 


38.5 


38. 1 


47.9 



NOTE. — Uncertainties on each spectral parameter are consistent with 
those reported for the tits in Table 3. 



these parameters are the temperatures (kT h ,kT w ; keV) and 
subtended angular sizes (f3h,f3 w ', degs) of the hot spot and 
warm ring, respectively, the viewing (ip) and hot spot (£) an- 
gles w.r.t the spin-axis (in degrees), the rotation phase (7; 
cycles), and finally, the NS radius (R; km) and distance (D; 
kpc). The XSPEC normalization is set to unity so that the flux 
is fixed by the distance and stellar radius, which implicitly 
takes into account all relativistic effects previously noted in 
§2. In the following spectral fits the pulsar distance is set to 
D = 3.3 kpc, based on radio pulse dispersion (Camilo et al. 
2006), and consistent with the measurement derived from HI 
absorption (Minter et al. 2007). 

This model allows us to predict the energy dependent mod- 
ulation, and use it to determine the viewing geometry and 
beaming pattern of the emitted radiation that best match the 
observations at different epochs. We assume both to be tem- 
poral invariant, i.e. no noticeable precession changes with 
time. In principle, the neutron star radius, because of general 
relativistic effects, could be uniquely determined; however, in 
practice the presence of noise does not allow for this most in- 



teresting of constraints. In the following, we apply this model 
to the data sets presented in Gotthelf & Halpern (2007), which 
fully described their preparation. We first present the spectral 
fits (§3.1) using the full model, assuming a face-on geometry 
(a = 0), and then constrain the overall emission geometry of 
the system (§3.2) by modeling the observed pulse modulation 
in 6 energy bands over time. 

All spectral fitting are done in the 0.7- 10 keV spectral band 
assuming no beaming initially, as this is not an important ef- 
fect spectrally. However, some degree of anisotropy of the 
radiation is found to be necessary in our model to reproduce 
the observed modulation (§3.2). 

3.1. Spectral analysis 

We started by fitting the phase-averaged XMM-Newton 
spectra for the 7 epochs simultaneously using our model for 
the pulsar emission geometry. Since the viewing geometry 
is not known a priori, we assume the simplest choice, that 
we are looking directly down the co-aligned rotation axis and 
magnetic pole (a = 0, see Eq. 1). This has the practical benefit 
of allowing the model code to run substantially faster since, 
for this special case of a = 0, the integration is simpler and 
only one call to the routine is needed for the computation of 
the phase-averaged spectrum. Across all epochs, all parame- 
ters are linked with the exception of the set of 4 epoch variable 
parameters (kT w ,kTh,f3 w ,(3h). Initial fits were used to deter- 
mine the nominal column density of Nh = 6.8 x 10 21 cm -2 , 
which was subsequently fixed to this value. 

An important technical issue for these fits is the degree of 
degeneracy between the radius R and the four epoch vari- 
able parameters (kT w ,kTh,p w ,Ph)- These 5 parameters over- 
determine the fit, unlike fits using a double blackbody model. 
Without fixing the radius there is no unique solution, and 
thus we consider a range of possible values between 9 < R < 
14 km, in 1 km increments. These results are presented in Ta- 
ble 1 and show a similar trend to those reported by Gotthelf 
& Halpern (2007) using the double-blackbody model. In both 
cases, the hot components is found to steadily decrease in size 
over time, while the warm component increases (with the ex- 
ception of the first data point). 

In our model, the radius of the star is not just a simple nor- 
malization. This is due to the introduction of gravitational 
redshift effects. Unlike for the non-relativistic case, the in- 
ferred temperatures of the spots increase as the radius of the 
star becomes smaller. Two counteracting effects, both due to 
flux conservation influence the spot size - gravitational red- 
shift tends to decrease the inferred emission area in the more 
relativistic (smaller) stars (due to the higher inferred tempera- 
tures); on the other side, for a fixed distance between the star 
and the observer, the spot angular size increases on smaller 
stars. For the values of the fit parameters here, the latter effect 
tends to dominate over the former. Over the sampled range, 
we do not find evidence for a preferred radius, based on the 
X 2 measurements. 

As discussed above, the results of our spectral fits using the 
above model show a similar trend to those reported by Got- 
thelf & Halpern (2007) using a double-blackbody model. The 
hot components steadily decrease in size over time, while the 
warm component increases (except for the first data point). 
We are aware of the importance of a possible third emission 
component from the rest of the NS surface, perhaps the qui- 
escent emission, initially masked by the significant extra flux 
from the warm component activated by the outburst. How- 
ever, we are unable to resolve any additional component, 
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TABLE 2 

Minimum xi as a Function of Beam Indices (n w ,nh) and NS Radius 



n w ,n h 






Minimum xt 








R=9km 


R=10km 


R=llkm S=12km 


ff=13km 


S=14km 


0,0 


4.06 


4.91 


3.97 4.62 


3.85 


2.86 


0.1 


0.86 


1.07 


1.13 0.96 


0.94 


0.99 


0.2 


3.17 


3.25 


3.74 3.02 


3.29 


2.83 


1.1 


2.51 


3.50 


2.55 2.64 


1.81 


1.17 



3.2. Modulation analysis 

Starting with the best spectral fit model parameters pre- 
sented in Table 1, obtained for £ = ip = 0, we now searched 
for the best values of £ and ip needed to reproduce the ob- 
served magnitude of the pulse modulation across the first 4 
epochs, as measured by Gotthelf & Halpern (2007). Given 
the uncertainties in the data, we limit our modeling to that of 
the pulse modulation, rather than the full pulse profile. The 
pulsed fractions were determined in six energy bands 5 at each 
epoch. For each value of the NS radius {R) fitted for in §3.1, 
we computed the modulation (defined below) over the grid of 
angles £,cf) < 90 deg, in 1 degree intervals. For each value 
on the grid, the model predictions were compared with the 
data. Notice that the flux depends on the angles £ and ip only 
through the parameter a in Eq.fHJ, and therefore it is symmet- 
ric with respect to an exchange of £ and ip. The magnitude of 
the model modulation is defined as 



PF = 



F +F 

1 max ' 1 m 



(13) 



In the geometry that we are considering, the maximum and 
minimum fluxes, F max and F m i n , correspond to phases 7 = 
and 7 = 7r, respectively. Both fluxes are integrated over the 
given energy bands. 

The results of these fits show that it is not possible to repro- 
duce the observed modulation if the emission pattern of both 
the hot and warm component is isotropic. The introduction 
of General Relativity effectively suppresses the modulation to 
below that observed, for any reasonable assumed NS radius. 
We ascribe the observed modulation to additional anisotropic 
emission from the thermal regions and test this assumption 
using a simple model of cosine beaming described in §2. For 
each value of the radius, we ran through the (ip,Q grid of mod- 
els for different combinations of the beaming parameters n w 
and «/, of the warm and hot components, respectively. More 
specifically, we varied «„,,«/, between and 2, in increments 
of 0.5 (note that the softer component is less modulated than 
the harder one). The fact that the warm component dominates 
in flux in the softest energy band, while the hot component 
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NOTE. — Minimum reduced xt a f ter comparing model and observed PFs 
over the £,»/>-space, with the NS radius and beaming indices held fixed at the 
given values. Only the first 4 epochs were included in this analysis. 



which is not required by the spectral fits. If the last few data 
sets are substantially affected by this potential third compo- 
nent, our pulse profile modeling of those data could be incom- 
plete. Therefore, for the second part of our analysis, we rely 
on the first 4 data sets alone, during which the emission from 
the two components dominates over that from the NS surface 
(whose quiescent level was measured with ROSAT prior to 
the outburst). 
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FIG. 2. — Reduced chi-square (xt) maps obtained by comparing modula- 
tion data and model described in the text for a range of viewing angles ip and 
£. The 68%, 90% and 99% confidence levels are shown for the best match to 
the observed pulsed fractions using the beaming patterns n w = 0, »j = 1, for 
R = 9 km and R = 14 km. The results are clearly degenerate with respect to 
an interchange of (£,ip). For the adopted model, the locus of minimum x'v 
depends only on radius and lies along a line (£ + /3(R))(ip + /3(R)) ~ constant 
(dashed lines). The minimum xi f° r eacn value of the radius is reported 
in Table 2, and the angles (ip,£) to which this minimum corresponds vary 
slightly with radius. 

dominates in the hardest energy band, allows us to constrain 
the degree of anisotropy of these two components indepen- 
dently. For each set of angles ip and £ on the grid, we com- 
puted the reduced xt an d k e Pt track of its minimum. Table 2 
reports the minimum xt tnat was obtained for a few represen- 
tative values of the beaming parameters. 

With the introduction of beaming, we are able to identify a 
set of model parameters that is able to reproduce the observed 
modulation in the first 4 epochs. The main results from our 
modeling can be summarized as follows: 

i) The modulation of the hot spectral component requires 
an anisotropic radiation pattern. For a cos"'' 5 emis- 
sion profile, the best match to the data is obtained with 
n h k 1. 

ii) No similar beaming is required to model the warm com- 
ponent modulation (i.e. n w ~ 0). 

Hi) We constrain the emission geometry by identifying 
allowed and forbidden regions in the ip-£ parameter 
space. 

iv) No NS radius is strongly preferred by the data. How- 
ever, the range of most preferred values of ip and £ 
varies with the radius of the star. 

For the optimal beaming parameters, n w = and «/, = 1, 
we plot the xt ma P computed for the analysis of Table 2. 
Figure |2] displays the 68%, 90% and 99% confidence lev- 
els drawn for the two extreme values of radii considered 
(R = 9, 14 km). As described above, this map is produced by 
comparing the model and observed modulation over a range 
of possible (£,ip) angle pairs, for our best fit spectral model 
parameters. The range of allowed solutions defines a locus of 
points in the £-?/;-plane along the line (£+P(R))(i/;+(3(R)) = c, 
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were c is, to zeroth order a constant for the given model (it 
varies by < 9% between the two extreme values of radii con- 
sidered here), while (3(R) contains the strong dependence on 
the NS radius. A linear regression fit yields the relations 
13 = 1.64i? km -6.16 and c = 40.65 R m + 1935, where R km is the 
radius of the star in km. The elongated shape of the contour 
plots shows that the two angles ip and £ are highly correlated 
in the fit. This is a result of the fact that the PF depends on a 
combination of these two angles. 

The dependence on the NS radius is clearly seen in Figure|2] 
For each given value of ip, the best fit loci moves toward larger 
values of £ as the radius of the star gets smaller. The contour 
levels for intervening values of the radii fall in-between those 
shown. This trend with radius can be understood as follows. 
As the radius decreases, the larger angular sizes and tempera- 
tures conspire to decrease the PF for the same viewing angles 
ip and £. In addition, the gravitational effect of light bending 
reduces the modulation even further for small stars. In order 
to reproduce the same observed PF for a given ip, a corre- 
spondingly larger £ is therefore needed for smaller stars since, 
for any given value of each of these two angles, a larger value 
of the other produces a larger modulation. This is the reason 
for the shift of the confidence levels toward larger values of £ 
for a given ip, when the radius of the star gets smaller. 

Figure [3] compares the data and model modulation for the 
(ip,0 values that yield minimum xt f° r tne case ^ = 12 km, 
riy, = 0, tif, = 1 . We show the results for two cases, one using 
all seven epoch data sets, and one for the case of the first four 
data sets only. In the former case, the minimum corresponds 
to ip = 49° and £ = 24°, while in the latter, it occurs for ip = 53° 
and £ = 23°. We find excellent agreement between data and 
model modulation for the first four observations alone; how- 
ever, the later data sets show increasing discrepancies, notice- 
ably increasing the overall contribution to \ 2 - This confirms 
the suspicion that, at later times, the emitted radiation acquires 
an unpulsed (or very mildly pulsed) contribution from the sur- 
face of the star. Our spectral fits, as described above, do not 
resolve this underlying stellar component, and therefore the 
spectral parameters close to quiescence might not be as repre- 
sentative of the underlying physical parameters of the system. 
However, the constraints on the viewing and emission geome- 
try that we derive using the first 4 data sets alone (cfr. Table 2 
and Fig. 2) can be considered robust, since the early data sets 
are basically unaffected by the presence of the star underlying 
emission. 

As a final step in our analysis we consider the validity of 
our initial method of assuming a face-on spectrum to derive 
the spectral model parameters that are then used to compute 
the modulation. As discussed in §3.1, the original fits were 
generated assuming ip = £ = for simplicity, prior to deter- 
mining the observational geometry of the pulsar system. We 
now show that this is an excellent assumption by refitting to 
the data the spectrum assuming the specific case of ip = 53° 
and £ = 23°, and then recomputing the modulations. The best 
fit spectral parameters are reported in Table 3 and shown in 
Figure|4] For these spectral parameters the model modulation 
is shown in Figure[3]as the dashed-line for the case of the first 
four data sets. The results are identical within the statistical 
uncertainty in the data. 

Similarly, we performed a test in order to assess the valid- 
ity of our method of analysis which separated the spectral and 
timing studies and used, as spectral parameters for the tim- 
ing analysis, those obtained from a phase-averaged fit to the 
spectrum. We again considered the specific case of ip = 53° 



and £ = 23°, and extracted phase-resolved spectra by dividing 
the observed spectrum into 5 equally spaced bins. We then 
fitted the model to the two bins centered on the maximum and 
minimum of the flux, respectively. We found that the temper- 
ature at flux minimum is lower than the temperature at flux 
maximum by < 10% for all XMM epochs with the exception 
of the first one, where the difference is ~ 20%. However, 
this difference is within one a of the combined uncertainty in 
temperature for all epochs. Therefore, our adopted method is 
quite robust within the statistical uncertainty of the data. 

4. DISCUSSION 

The time-dependent spectrum and pulse modulation of the 
transient magnetar XTEJ18 10-197 provide a unique diag- 
nostics of its emission properties and geometry. Under the 
assumption that the post-burst emission is described by two 
thermal components, as early analysis of this object sug- 
gested, we have been able to extract information on some of 
the physical properties of the star, through a detailed modeling 
of the combined spectra and pulsed modulation together. 

We found that, while the phase-averaged spectral fits alone 
are degenerate with respect to the emission pattern of the radi- 
ation, including modeling of the energy-dependentpulsed flux 
allows us to constrain the properties of the emission region. In 
particular, since the warm component dominates in the lowest 
energy band, while the hot component dominates in the high- 
est energy band, the PFs are able to determine the degree of 
anisotropy of these two components independently. We found 
that the warm component is best described by an isotropic 
emission pattern, while the hot component is well represented 
by an emission pattern of pencil type, f(8) oc cos S, where S 
is the angle that the emitted photons make with respect to the 
normal to the surface of the star. The different type of radi- 
ation pattern required by the low and the high energy com- 
ponents could be seen as an indirect confirmation of our as- 
sumption that the contribution from these two energy bands 
does indeed come from different components. Beaming of ra- 
diation in the direction of the magnetic field is predicted by 
models of magnetized atmospheres in the limit of high mag- 
netic fields (e.g. van Adelsberg & Lai 2006). Since the hot 
spot is produced in a region much smaller than that of the 
warm component, it is more likely to find a configuration with 
parallel field lines in the hot region (and most likely perpen- 
dicular to the surface, which favors the heat flow), than in the 
warm region. The latter might rather encompass regions with 
different orientations of the magnetic field, hence resulting in 
an overall more homogeneous radiation pattern. 

The strongest constraints that we derived from our analysis 
are on the geometry of the star. In particular, we determined 
the allowed regions for the angle (£) between the spot axis 
and the rotation axis, and the angle (ip) between our line of 
sight and the rotation axis. These two angles determine the 
minimum and maximum angles between the line of sight and 
the spot axis, given respectively by a m \ n = £ - ip and a max = 
£+ip. We find that, while the range of a m \ n is compatible with 
very small angles (including zero), however a max must always 
be large, > 60° within 3a confidence level for any value of the 
star radius. Being able to rule out to a high confidence level 
small viewing angles a for the entire rotation period of the 
star bears important implications for models of the observed 
radio emission from this object. In fact, Camilo et al. (2007a) 
showed that the peaks of the radio and the X-ray pulses are 
aligned, suggesting that the footpoints of the active magnetic 
field lines on which radio emission is generated are also the 
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FIG. 3. — Comparison between the measured pulsed fraction (PF) and that predicted by our model for one particular set of best-fit parameters determined for 
a R = 12 km neutron star. Results are shown for the best fit model using all seven data sets (red line; ip = 49° and § = 24°) and for the first four data sets only 
(green line; ip = 53° and £ = 23°). At the later times, the model is seen to deviate significantly from the data. A likely explanation is the increased contribution 
of the unmodeled emission from the stellar surface over time compared to the modeled flux. In the first 4 data sets, the dashed line is to be compared to the green 
line, showing the effect of using the face-on spectra instead of the iterated spectra for the given viewing geometry (see text for details). Each data point is drawn 
in the middle of the corresponding energy band, for those bins with sufficient photons. 



TABLE 3 

Spectral Results for Radius R = 12 km and Viewing Angles ip = 53 deg, and § = 23 deg; D = 3.3 kpc, Nh = 6.8 x 10 21 cm 2 



Parameter 


Sep 8 


2003 

Oct 12 


Mar 11 


2004 

Sep 18 


Mar 18 


2005 

Sep 20 


2006 
Mar 12 


kT h (keV) 
kT„ (keV) 
Ph (deg) 
Pw (deg) 


39+| 


a87 S?t 

34+" '' 
'■ L -Q ! 2 

30+2 


0.85+00' 

J - J -o.i 
31.0^ 


a82 1i 

4 5j 2 +^ : ^ 
^■ JZ -0.08 

33.9+0;^ 


73+0-03 
-002 
26+0 5 

T. S+0.4 
3 - S -0.2 

37 7+°- 7 

* —0.1 


63+0 03 
-Q.04 

23+0 0' 
-Q.03 

J - uz -0.05 

44+' 


53+0- 04 
-Q.Q3 

()21+oo 2 
60+5 



NOTE. — Uncertainties in spectral parameters are 90% confidence for two interesting parameters. The x 2 of the fit is 1.08 for 1692 dof. 



locations of the concentrated heating that is responsible for 
the enhanced X-ray emission, This means that, even if the 
radio emission is likely produced at much higher altitudes on 
the surface of the star than the X-ray emission, however the 
axis where the two emissions peak is the same (or very close). 

Attempts to constrain the viewing angles of 
XTEJ18 10-197 using radio polarimetry were made by 
Camilo et al. (2007b). They found two configurations 
likely, one with £ ~ 70° and a m [ n ~ 20° -25°, and another 
with £ ~ 4° and a m ; n ~ 4°. Our fits rule out the second 
configuration to a high significance level. The value of 
ctmin ~ 20° -25° on the other hand is perfectly compatible 
with our results, albeit it requires £ ~ 60° if R = 9 km and 
£ ~ 50° if R = 14 km. Although our confidence levels are 
close to one of the two solutions of Camilo et al. (2007b), 
we cannot make a formal statistical comparison with their 
results, since they do not have a reliable estimate of the 
parameter uncertainties from their radio measurements (F. 



Camilo, priv. comm.). 

As discussed by Camilo et al. (2007b), the observed wide 
radio pulse profile of ss 0.15 P can be explained by either a 
model in which the magnetic and rotation axes are almost 
aligned, or by a model in which the emission height is very 
large. Our results strongly rule out the first scenario, hence 
implying a large emission height. This, in turn, implies a 
large opening angle of the beam (Gil et al. 1984), compa- 
rable to that observed in young pulsars (Johnston & Weisberg 
2006). These characteristics of the radio emission, if common 
in magnetars, make more stringent the limits on the radio for 
the greatest majority of the objects that have not been detected 
in this waveband, and leave even more open the question of 
what is that makes some magnetars different. 



We thank Jules Halpern, Andrei Beloborodov and Fernando 
Camilo for stimulating discussions on several aspects of this 
work. RP thanks Columbia University for the kind hospitality 
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FIG. 4. — XMM-Newton phase-averaged spectra of XTE J 1 8 1 0— 197 obtained at 7 epochs fitted with the model presented in the text, for the specific case of 
R = 12 km, tj) = 53° and £ = 23°. The parameters for the best fit model are reported in Table 3. The lower panel shows the collected residuals to this fit for each 
spectrum. 

during the several visits made while this work was carried out. comments on our manuscript. 
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